# NOT RUN {
## 1 - GEV fit
x <- rep(NA, 100)
for (i in 1:100)
x[i] <- max(rnorm(365))
gevmle(x)
## 2- GPD fit
x <- rnorm(10000)
##we need to fix a threshold
u <- quantile(x, 0.99)
gpdmle(x, u)
# }
Run the code above in your browser using DataLab